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Abstract 

A 3-parameter family of helical tubular surfaces obtained by screw 
revolving a circle provides a useful pedagogical example of how to study 
geodesies on a surface that admits a 1-parameter symmetry group, but is 
not as simple as a surface of revolution like the torus which it contains 
as a special case. It serves as a simple example of helically symmetric 
surfaces which are the generalizations of surfaces of revolution in which 
an initial plane curve is screw-revolved around an axis in its plane. The 
physics description of geodesic motion on these surfaces requires a slightly 
more involved effective potential approach than the torus case due to 
the nonorthogonal coordinate grid necessary to describe this problem. 
Amazingly this discussion allows one to very nicely describe the geodesies 
of the surface of the more complicated ridged cavatappi pasta. 
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Figure 1: The ridged cavatappo (singular) / cavatappi (plural) pasta shape. 

1 Introduction 

In looking for examples of interesting surfaces to use in playing with geodesies, 
one starts with the plane and cylinder (flat, rotationally and translation sym- 
metric), moves on to the sphere (rotationally symmetric), whose great circle 
geometry was long used in celestial navigation even before non-Euclidean geom- 
etry and then differential geometry found them to be nice mathematical objects 
of study. These are all maximally symmetric (both isotropic at every point: 
every direction is equivalent, and homogeneous: each point is equivalent to ev- 
ery other point) with 3-dimensional symmetry groups. After this, surfaces of 
revolution (which have rotational symmetry about an axis) present themselves 
as candidates for study, with the 2-parameter family of tori as the first inho- 
mogeneous but still rotationally symmetric simple surface one encounters next 
yet still constructed from circles, demonstrating both positive and negative cur- 
vature, and containing the sphere as a special case. This provides a nice test 
case to study geodesies using the physics approach to motion along curves con- 
strained to lie on a surface, revealing lots of fun mathematics developed in detail 
in a previous investigation pQ. 

What next? It turns out that the tubular surface based on a helix is a 3- 
parameter family which contains the torus as a special case where the helix de- 
generates to a circle, but it still has a 1-parameter symmetry group of corkscrew 
rotations (a linked rotation about and translation along the central axis of the 
helix) that makes its study reducible to motion in 1-dimension, like the torus. 
In fact amazingly many of the properties of geodesic motion on the torus apply 
to the helical tube. One can imagine constructing one screw-revolution of this 
tube by cutting a horizontal torus by a half-plane through its vertical axis of 
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symmetry along a vertical circular cross-section (a meridian) , and then stretch- 
ing the two ends vertically in opposite directions until the circle of centers of all 
such cross-sections stretches into a segment of a helix about the axis of symme- 
try. One can then extend the terminology of inner and outer (helical) equators, 
(circular) meridians, (helical) parallels, and northern and southern polar helices 
in a natural way from the corresponding curves on the torus, in turn taken 
from the sphere. Remarkably the equators remain geodesies of the surface as 
in the torus case, while the study of general geodesies can be reduced to the 
1-dimensional problem of "radial" motion along the meridians due to the con- 
served "screw-angular momentum" associated with the 1-parameter symmetry 
group, generalizing the conserved vertical component of angular momentum of 
the corresponding torus case. The terminology of bound and unbound geodesies 
also generalizes to the new case, with a qualitative correspondence between the 
bound geodesies oscillating about the outer equator but not crossing the inner 
equator, and those unbound geodesies which do cross the inner equator. 

This larger class of surfaces containing the tori is useful in that it adds a 
new element to the discussion: a nonorthogonal natural coordinate grid on the 
surface. This is easily orthogonalized in 2-dimensions by completing the square 
to adapt a new basis of the tangent spaces to the orthogonal decomposition 
with respect to the parallels which are the helical orbits of the 1-dimensional 
symmetry group of this geometry. Once done, all of the qualitative features 
of geodesic motion on tori generalize to the larger family of helical tubular 
orbits. We are all familiar with this symmetry: turn a metal screw in its hole 
by 360 degrees and it returns to its original orientation except moved farther 
into the hole. The angle of rotation is proportional to a translation along the 
axis of the rotation. The present discussion is understood to be a sequel to the 
torus article pQ , showing how to take into account the nonorthogonal coordinate 
system to reduce the geodesic discussion to one qualitatively like the simpler 
torus case. 

For particular values of the shape parameters of our family of helical tubular 
surfaces, two and a half turns of the tube models extremely well the ridgeless 
(liscio) cavatappi pasta shape (also known as gobetti, spirali or cellentani), 
which is a common corkscrew version of the macaroni family of tube pastas 
whose design aspects were recently documented by George Legendre in this 
elegant book Pasta by Design [2j[3]. The default "rigato" (ridged) cavatappo 
that is actually found in markets is a slight generalization of the helical tubular 
shape within the class of surfaces admitting a helical symmetry and although 
its geodesic equations require several pages to write out, they are easily handled 
numerically by the computer algebra system required to derive them. For this 
pasta shape, this exercise takes the mathematics of pasta one step beyond the 
design aspects documented by the architect George Legendre to the dynamics 
of its intrinsic geometry as approached by a physicist. 

Although the screw-rotation symmetry is familiar to most of us through its 
everyday representation in bolts and metal sheet screws (but not wood screws 
which taper off) and cork screw wine bottle openers, the corresponding mathe- 
matical symmetry group is surprisingly absent as an example in mathematical 
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Figure 2: An illustration of the inclination angle 77 from the horizontal of the 
central helix of the surface in <f>-z plane of the tangent space at a point along 
that curve, whose unit tangent vector is aligned with the hypotenuse of the right 
right triangle, and binormal with the hypotenuse of the left right triangle, while 
the normal — p points along the cylindrical radial coordinate line towards the 
symmetry axis, where (/5, 0, z) are unit vectors along the cylindrical coordinate 
lines. The value of the inclination angle 77 = arctan(8/15) « 28.1° illustrated 
here is very nearly equal to Legendre's choice for his cavatappo design parameter 
but corresponds to the fourth Pythagorean triple (8, 15, 17) right triangle leading 
to rational values of the surface parameters: (a, b, c) = (|, 1, |) = | (l, |, ^) . 

discussions of continuous transformation groups [4] or in toy applications in clas- 
sical mechanical systems. However, the Archimedes screw has been with us for 
millcnia and still finds useful applications in large scale industrial screws that 
benefit society in various ways [5]. Surprisingly, a careful mathematical study 
of improving the traditional design was only done recently by Chris Rorres with 
the help of a computer algebra system [H| , with a followup study resulting in an 
even more optimal design. 

2 Helical tubular surfaces 

Starting from the usual Cartesian coordinates (a;, y, z) in Euclidean space where 
the element of differential arclength ds satisfies an iteration of the Pythagorean 
theorem ds 2 = dx 2 + dy 2 + dz 2 (the "line element" representation of the "met- 
ric" of Euclidean space), a simple tubular surface built from a helix can be 
represented by the following parametrization [2l|3] 

x = (a + b cos x) cos </> , y = (a + b cos x) sin cf) , z = c (f> + b sin \ > (1) 
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or equivalently as p = a + b cosx, z — c<f> + b sinx i n the usual cylindrical 
coordinates (p, <fi, z) related to the Cartesian ones by x = p cos <f>,y = p cos <f>,z = 
z. It is assumed that a > b > 0, and for convenience, that c > (for a right- 
handed tube, c < for a left-handed tube). The surface lies between the inner 
and outer cylinders of radius a — b<p<a + b around the symmetry axis taken 
here to be the vertical z-axis. Setting c = reduces this surface to a torus. 
By restricting the two angles x an d <j> to either interval (— ir,w] or [0, 2w) as 
convenient, they may be interpreted as coordinate functions on the surface. 

The "central" helical curve corresponds to setting the radius of the verti- 
cal circular cross-section of the tube to zero: 6 = 0, while a is the radius of 
the cylinder containing that central helix, and c is its inclination parameter, 
with "coiling" angle of inclination r) = arctan(c/a). The geometry of this curve 
parametrized by <fi is easily seen by introducing the unit vectors along the cylin- 
drical coordinate lines 

p = (cos (j>, sin <j>, 0) , 4> — {— sin0, cos</>, 0) , z = (0,0,1), (2) 

in terms of which it and its tangent have the following representation 

r{<f>) = ap + c(f)z , r '((f)) = acf) + cz (3) 

illustrated in Fig. [5J 

Fig. [3] illustrates the construction and one complete revolution of a helical 
tubular surface with its inner (\ = vr) and outer (\ — 0) equators marked 
off. The grid shown in the computer rendition of the surface consists of the 
constant 4> circles which result from the intersection of the torus with vertical 
planes through the symmetry axis (the meridians: cj> = </>o) and the constant 
X = Xo helices (the "parallels," similar to the parallels in surfaces of revolution). 
The northern (x = 7r /2) and southern (x = —n/2) polar helices correspond to 
the northern and southern polar circles on the torus which in turn generalize the 
north and south poles on the sphere. The "radial" arc length coordinate r = bx 
and the corresponding angle x are measured upwards along the meridians from 
the outer equator. One can also introduce the four hemitubes generalizing the 
corresponding hemitori of the torus: the inner (—tt/2 < x < ' r /2), the outer 
(tt/2 < x < 37t/2), the upper (0 < x < n ), an d the lower (— n < x < 0) 
hemitubes. 

If one thinks about how to generate a tubular helical surface with circu- 
lar cross-sections, there is no reason why not to consider a circular cross-section 
that is tilted with respect to the vertical direction assumed in the present model 
(cavatappo 1.0). It is easy to add one additional parameter to this family of 
basic surfaces which is the angle of tilt of the normal to the plane of the circu- 
lar cross-section in the 4>-z plane of the tangent space up from the horizontal 
direction. In this enlarged family (cavatappo 2.0), there is an obvious preferred 
tilt angle which makes the plane of the circular cross-section coincide with the 
normal plane to the central helical curve. This new orthogonally tilted cavat- 
appo model has very nice mathematical properties which are discussed in a 
subsequent article [7]- 
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Figure 3: Above: A vertical half-plane cross-sectional circle of the helical 
tubular surface built around a helix through the center of this circle whose axis 
of symmetry is the z-axis. This circle in the x-z plane is simultaneously rotated 
around this axis while being translated upwards along that axis (c > 0) , so that 
the right hand rule wrapping fingers around the helix in the direction in which 
it is rising (lower figure) puts the thumb up. Below: This figure illustrates one 
turn < 4> < 2ir of the "unit tube" case (a, b) = (2, 1 ) of a unit circle which 
is revolved and translated around and along the z-axis with inclination angle 
?7 = arctan(2/5) « 21.8°, with an inner equator always a unit distance from the 
axis, a cavatappo surface model due to Sander Huisman [3]. The outer (\ — 0, 
red) and inner (\ = ^ zn J green) equators are shown together with the "prime 
meridian" (9 — 0, blue). The northern and southern polar helices (x — ±tt/2, 
cyan) are equidistant from the two equators. 
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Substituting the differentials of these coordinates into the Euclidean metric 
ds 2 = dx 2 + dy 2 + dz 2 to evaluate the induced metric on the surface, one finds 
easily for its line element representation 

ds 2 = ((a + 6cosx) 2 + (?) dcj) 2 + 2cb cosx d<f)dx + b 2 dx 2 ■ (4) 

This metric is independent of <f>, which means that it is invariant under trans- 
lations of the azimuthal angle <f> along the family of parallels, termed helical 
symmetry [4] or more suggestively "screw-rotation" symmetry. The vector field 
d/d(f> on the surface which generates these translations is said to be a Killing 
vector field of this metric. For motion along geodesic curves within the surface, 
the component of its affinely parametrized geodesic tangent along the Killing 
vector field remains constant, a "contant of the motion," or a conserved mo- 
mentum associated with this symmetry group. This Killing vector field is just 
the restriction to the surface of the Killing vector field of the Euclidean metric 
Cspacc = xd/dy — yd/dx + cd/dz — d/d<f) + cd/dz expressed in either Cartesian 
or cylindrical coordinates, which generates a rotation about the "screw axis" by 
an amount which is proportional to a simultaneous translation along the di- 
rection of that axis. In cylindrical coordinates these corkscrew rotations are 
p — > p, 4> —> <p + t, z — » z + ct. Within the surface, these are just translations in 
4>, so the Killing vector field expressed in the surface coordinates, namely the 
intrinsic Killing field, is just £ = d/d(j>. 

If we introduce the arclength radial coordinate r and cylindrical radius func- 
tion R(r) for the helical center curve 

r = b\ 1 R = a + bcos(r/b) , (5) 

the metric can be written in the following form as well as a second form obtained 
by completing the square on the differential d<j>, for which we give two versions in 
order to highlight the conserved momentum combination we shall see emerging 
for geodesic motion 

ds 2 = (R 2 +c 2 )d(f> 2 + 2ccos(r/b)d<pdr + dr 2 

,_ 2 a ./,, ccos(r/b) \ 2 ( R 2 + c 2 sin 2 (r/6) 

((R 2 + c 2 )d4> + c cos(r/&) drf / R 2 + c 2 sin 2 (r/fr) \ 
R 2 + c 2 + \ R 2 + c 2 ) 

Completing the square adapts the metric to the orthogonal decomposition of 
the tangent space with respect to the intrinsic Killing vector field £ = d/d(f). 
This metric is invariant under reflections through the origin of the coordinates: 
(r,tj>) -> (-r, -<f>). 

Note that the metric determinant (from the second line, an orthogonal frame) 
is the factor R 2 +c 2 sin 2 (r/6), whose square root provides the surface area density 

1 /2 

function: dS = (i? 2 + c 2 sin 2 (r/6)) dr d4>- This can be integrated exactly over 
one revolution of the surface, but the result is an extremely long formula infested 
with elliptic functions which has little value in being reproduced. 



J dr 2 (6) 
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If (r(A), </>(A)) is an affinely parametrized geodesic of this metric on the 
surface, with tangent 

u = i\ dr + = ur£>r + v * d + ' (7) 

then its orthogonal decomposition is 

U = U f e? + U^e^ (8) 
expressed in terms of its components 



U r = 



R 2 + c 2 sin 2 (r/6)\ 1/2 dr 



R 2 + c 2 J dX 



with respect to the orthonormal frame 



R 2 + c 2 sin 2 (r/b)\ 1/2 ( ccos(r/6) 



R 2 + c 2 J V i? 2 + c 2 



whose dual frame is 

?2 1 „2 c ,:„2/„ \ 1 / 2 
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(ii) 

The component of the tangent vector along the Killing vector field is a 
conserved screw- angular momentum, i.e., a constant along the geodesic 

> = w u = (R2+ ^{t) +c «w»> (I) ^ < fl2 + < 2 > 1/2 °* • < 12 > 

which is exactly the combination occurring in the last form of the metric adapted 
to the orthogonal decomposition with respect to the radial direction along the 
meridians. In an affine parametrization, the square of the length of the tangent 
is also a constant, which we will call twice the energy 

dsV t 2 / fl 2 +cW(r/&) \ (drV 

dx) = WT^ + { — WT7 2 ) [dx) =2E ' (13) 

where we have substituted for d(f>/dX in terms of the screw-angular momentum 
£. The length of the tangent vector in an affine parametrization of a geodesic is 
a constant along that geodesic, with E = | for an arclength parametrization in 
which this length is 1. 
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If we interpret this problem as geodesic motion in the surface in the physics 
language of motion in space where the affine parameter A plays the role of the 
time (and U is then called the velocity vector and U r ,U^ the velocities), E 
and i are called "constants of the motion." Since both £ and E are constants of 
the motion, we obtain a single constraint on the square of the "radial velocity" 
dr/dX, or equivalently the ortho normal component U r , which is of the form 

±(U f ) 2 + V = E, (14) 

where 

V ~ 2(R 2 + c 2 ) ~ 2 ((a + bcos(r/b)) 2 + c 2 ) (15) 

acts as an "effective potential" for the radial motion alone. Note that it is 
an even function of the radial variable r, as well as periodic with period 2irb. 
Qualitatively, this potential leads to the same kind of radial motion as for the 
special case c = of a ring torus thoroughly discussed in the previous article 
PP. One can fix the potential for i ^ by choosing £ to have a particular 
value, and then using the freedom of the affine parametrization to vary the 
energy E corresponding to energy levels in the graph of the effective potential 
to describe the allowed interval of radial motion where V < E. When V = E at 
r = r± corresponding to dr/dX = 0, turning points of the motion occur. These 
intervals r_ < r < r + to which such motion is confined are symmetric about 
r = 0: r_ = — r+. The outer equator and the corresponding geodesies may 
be called bound orbits in analogy with the torus problem and more generally 
1-dimensional motion in a potential well. If E > V, the motion is unbound with 
unbounded values of the radial coordinate corresponding to an infinite number 
of crossings of the inner equator. Local minima of the potential then correspond 
to stable equilibria, while local maxima correspond to unstable equilibria. The 
outer and inner equators are themselves geodesies which correspond respectively 
to a stable and unstable equilibrium. 

Note that solving the energy equation for dr/dX = f(r), the problem is 
reduced to a quadrature 

A= f f(t)dt (16) 

which must be evaluated numerically, after which the screw-momentum equation 
can be rewritten for (f> = 4>(r) 

I = F(r, dr/dX, d<f>/dX) = F{r, f(r)-\f(r)d<j>/dr) (17) 

and solved for dip/dr and integrated formally to give <p as a function of r (only 
amenable to numerical integration) , but this gives no overview of the classifica- 
tion of the orbits which result from the process. The effective potential instead 
gives a simple visual description of that classification and enables particularly 
interesting orbits to be studied. 
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The second order geodesic equations are 




Notice that (dr/dX,d<p/dX) — > (~dr/dX, —d<fi/dX) is a symmetry of this system 
which maps £ — > —£, and so its initial data at a given initial position is reflection 
symmetric about the origin in the initial tangent space. 

Note that r — 0, bir reduces the geodesic equations to d 2 <f)/dX = requiring 
the azimuthal angle to be linear in the afhne parameter A along the inner and 
outer equators, which are therefore geodesies as in the torus case. 

If one evaluates the Gaussian curvature function on the surface (half the 
Riemann scalar curvature) using a computer algebra system, one finds 

1 cos(r/6)(a + 6cos(r/6)) 3 — be 2 sin 4 (r/&) , , 

K = 7j . (19) 

b ((a + 6cos(r/6)) 2 + c 2 sin 2 (r/6)) 

The sign of the curvature would be the sign of cos(r/6) if it weren't for the ex- 
plicitly negative term in the numerator, which pushes the zero of curvature from 
the polar helices towards the outer equator, surrounded by a zone of positive 
curvature like the torus. 



3 Generalized helical surfaces: the real cavat- 
appi surface as an instructive example 

One may generalize this family of helical tube surfaces by letting the circular 
meridian be replaced by any closed plane curve not intersected by the screw-axis, 
but one needs to be a little more organized about the corresponding calcula- 
tions. This is essentially an application of the threading splitting approach to 
a general relativistic spacetime and so serves as a useful example for practice 
with the associated ideas in a more concrete setting [5] . An instructive example 
is the actual cavatappi pasta shape, which is "rigati" (ridged) by a small sinu- 
soidal perturbation of the circular cross-section, with amplitude d and integer 



10 




Figure 4: A vertical half-plane cross-sectional circle of the cavatappi (rigati) 
pasta shape, together with one turn of the surface. The radial coordinate r = b\ 
is no longer an arclength coordinate on the surface, but instead measures the 
arclength along the original circular cross-section. A bound geodesic is shown 
which makes a bit more than 11/2 oscillations about the outer equator. 
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frequency m > 1 (an integer is necessary for periodicity in turn required for 
the surface to be smooth), leading to to — 1 outward bulges and to — 1 inward 
bulges. This is a 5-parameter family of surfaces described by the equations 

x = (a + bcos\ + dcosmx) cos <j> , 

y = (a + bcos\ + dcosmx) sin (j) , (20) 

z = ccfi + bs'mx + dsinTO-x , 

with — oo < (f), oo, — <j) < x < <l>- The pasta design expert George Legendre [2] 
gives the parameter values (only the ratios of the first four are important for 
the shape, we scale his values down by a factor of 2): 

a = l, b = l, c=l,d=±r, m=10, (21) 

2 5 21) 

one complete revolution of which is shown together with its meridian cross- 
section in Fig. |U The pasta shape has 2.5 complete revolutions: < < 5ir. 
Marked also is a bound geodesic path starting at the outer equator. Actually 
George uses c = 5/2it k 0.796 « 4/5 = 0.800, but the nearby rational value of 
the parameter slightly simplifies formulas. 

If one looks closely at the cross-section in Fig. |4l one sees an asymmetry 
between the outward bulges and the inward depressions. A symmetrical ridging 
would have first set d = and instead replaced b by b + dcos(mx) in Eq. (fTJ). 
Clearly one can study many variations of this theme with the same techniques 
used below, nor do the particular choice of cavatappo parameter values have 
any more than an approximate meaning in the real pasta world. 

If in addition to r = bx we introduce a few abbreviations 

R = a + bcos(r/b) + dcos(mr/b) , 

S = bcos(r/b) + dmcos(mr/b) , (22) 

T = 6sin(r/6) + dmsin(mr/b) , 

the metric can be written in the following orthogonal form 

ds 2 - (R 2 + c 2 ) dcp 2 + 2^ # dr + ( ^ ± r2 ) dr 2 

= <?00 dcf) 2 + 2g<f, r d(p dr + g rr dr 2 (23) 
= M 2 (d(p + M r dr) 2 + ^ rr dr 2 , 

= (a/) 2 + (u/ ) 2 , (24) 
which requires computation of the following quantities: 



M = (g M ^ 2 = (R 2 + c 2 Y /Z , M, 



1/2ll5 5 2 + T 2 -c 2 5 2 /(i? 2 + c 2 ) 

7rr = grr ~ M M 2 = — ^ . (25) 



<?0 6(i? 2 + c 2 ) ' 

v 

6 2 
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This corresponds to the orthogonal frame and dual frame 



e r = d r - M r d<f>, e<p = d lt> 

uj r =dr, uj^ = d<t> + M r dr (26) 
and orthonormal frame and dual frame 

ef = y rr / e r , = M I 

u r = 7r 1 / 2 w r , ^ = M^ . (27) 
Note that MjU 2 = {R 2 + c 2 ){S 2 + T 2 ) - c 2 S 2 is the determinant of the metric 

1 /2 

and determines the surface volume element M^ rr dr d<p of the tubular surface. 

If (r(A), 0(A)) is an affinely parametrized geodesic of this metric on the 
surface, its tangent is 

U = ^d r + ^ = lfe f + U*e $ (28) 
and has orthonormal components 

These can be inverted to yield 

U r = ^-= jrr^U' , U* = ^= M~ l U* - M rl -V 2 U f . (30) 
dX dX 

The component of the tangent vector along the Killing vector field (the 
conserved screw-angular momentum) is a constant along the geodesic 



£ — MU^ — M 2 (%+Mr£) 
\dX dA J 



(31) 



as is its length, half of which we call the energy 



so again 



1 I 2 1 I 2 
V = ___ = _ (33) 

2M 2 2 (a + bcos(r/b) + dcos(mr/b)) 2 + c 2 

acts as an effective potential for the radial motion. If the screw-angular mo- 
mentum is nonzero, then we might as well set it equal to 1 and use the energy 
parameter to distinguish initial data. Fig. [5] shows the nonzero screw-angular 
momentum potential for the cavatappi pasta surface together with the corre- 
sponding smooth shape with d = 0. 

When I ^ 0, of course we can solve the energy constraint on the radial 
velocity for that derivative dr/dX = f(r) and then integrate dX = dr/f(r), but 
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Figure 5: The orthogonal decomposition of the radial direction with respect 
to the parallels in the tangent space (top). The initial data is specified as a 
unit vector with angle (3 measured down from the direction orthogonal to the 
parallels (bottom). 
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Figure 6: The effective potential for geodesic motion on the cavatappi pasta 
surface (solid) and its smoothed version (dashed). Since the frequency m = 10 is 
even, there are an odd number of bulges, so in addition to the the inner equator 
geodesic, there is a pair of unstable geodesies on either side of it on the two 
nearby innermost bulges closest to the screw axis. 



this is extremely complicated and must be done numerically. If i = 0, then we 
get a simple first order differential equation to integrate for the angle as a 
function of the radial variable r 

^ + M r $- = 0^d<f>= -M r dr -> 
d\ d\ 

f r c 

M r (w)dw = — / - (bcos(w/b) + dm cos((mw / b)) dw . (34) 
Jo " 

In particular one loop of the surface corresponds to 

A<f> = - M r {w)dw ^0. (35) 
Jo 

which does not vanish since this is a function which is even about the inner 
equator <f> = n which is more negative near the inner equator than it is near the 
outer equator <j> — 0. Thus these are not closed geodesies. 

If we introduce polar coordinates in the tangent plane ("velocity space") 
adapted to the orthonormal frame, then the velocity has constant length {2E) 1 / 2 
because of the affine parametrization assumed for the geodesic, so we get the 
representation 

(U f ,U^) = (2S) 1 / 2 (cos/?, sin/?) (36) 
in terms of the angle measured from the direction orthogonal to the parallels, 
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which in turn implies the constancy of 

^A^Msin/3. (37) 

For a given energy that sets the afhne scale, and for a given initial position on the 
surface, the initial angle /3 parametrizes the family of geodesies which emanate 
from this point. If we choose to numerically solve the coordinate component 
form of the geodesic equations rather than the orthonormal component form, 
then Eq. (|30[) expressed in terms of /3 gives us geometric initial velocity data 
for those second order differential equations. Setting E = 1/2 then corresponds 
to an arclength parametrization. If instead we set £ = 1, then the relation (|37[) 
sets the energy level according to the initial values of r and f3. 

At this point one can study the details of this system in the same approach 
taken in the preceding article [1] written to be understood by undergraduates 
who have already taken multivariable calculus, differential equations and linear 
algebra. A computer algebra system is essential to make concrete calculations. 
While George Legendre chose Mathematica, the author is a long time Maple 
user. 



4 Italy and Tensor Calculus/Differential Geom- 
etry 

Dirk Stuik also recalled that when asked what he liked best about Italy, Einstein 
responded "spaghetti and Levi-Civita." 
—quote from a talk at a 1996 AMS Meeting [9]. 

While there is widespread awareness of the association between pasta and 
Italy, few people realize the key role Italian mathematicians played in develop- 
ing the tools of tensor analysis and differential geometry which made Einstein's 
theory of general relativity possible. While the initial ideas of curvature of 
surfaces were developed by Gauss, and Riemann generalized them to arbitrary 
dimension, it was Gregorio Ricci-Curbastro who provided the mathematical 
infrastructure of tensor calculus to make it usable for calculating, and his stu- 
dent/collaborator Tullio Levi-Civita [KJ] joined this enterprise, codifying the 
ideas into a landmark article in 1900 known for a time as Ricci Calculus [TTj . 
Einstein was not very mathematical, and it was this article that his friend Mar- 
cel Grossmann [T2] used in 1912-1914 to teach him enough to get serious about 
Einstein's attempts to generalize special relativity, and Levi-Civita participated 
in Einstein's progress through and after the birth of the final theory of general 
relativity published in 1916, studying its solutions and introducing the key idea 
of parallel transport in 1917 [13]. This latter concept plays a fundamental role in 
both using general relativity and in understanding differential geometry. Levi- 
Civita presented this theory in lecture notes put together with Enrico Persico 
(friend of Enrico Fermi, and editor of his collected works) in 1925 and translated 
into English in 1927 as The Absolute Differential Calculus [T4" l ll5 j . 
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The present family of surfaces provide an interesting playground for students 
to learn a bit of geodesic theory while perhaps having some fun. Yet unexplored 
is the extrinsic curvature properties of these surfaces. However, the discussion 
so far also shows how essential computer algebra systems are for actually doing 
the calculations once some intelligent framework is imposed on the problem. 
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